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Abstract. The conformational properties of flexible polymers in d di- 
mensions in environments with extended defects are analyzed both 
analytically and numerically. We consider the case, when structural 
defects are correlated in Sd dimensions and randomly distributed in 
the remaining d — Ed- Within the lattice model of self-avoiding ran- 
Ch ' dom walks (SAW), we apply the pruned- enriched Rosenbluth method 

Ch , (PERM) and find the estimates for scaling exponents and universal 

' shape parameters of polymers in environment with parallel rod-like 

defects {sd = !)• An analytical description of the model is developed 
within the des Cloizeaux direct polymer renormalization scheme. 

» 

1 Introduction 

Ch ■ Conformational properties of polymer macromolecules in solutions are the subject of 

^ I great interest in statistical polymer physics |ll2l3j . It is established, that long flexible 

polymers are characterized by a number of universal characteristics, independent 
of the chemical structure and architecture of the molecules. A typical example of 
such a quantity is the averaged end-to-end distance of A^-monomer polymer chain 
' (characterizing the linear size measure of macromolecule), that scales according to: 

o 



(N 



X 



(R^) - N^", (1) 



where z/ is a universal exponent, depending on space dimension d only {v{d = 3) = 
0.5877 ± 0.0006 Q). An other example is the rotationally invariant characteristic of 
the shape of typical polymer conformations, namely the aspherisity (^4^) [S], which 
equals for sphere-like shapes and 1 for rod-like. For a polymer in a pure solvent in 
d = 3 it was found (^3) = 0.431 ± 0.002 0. 

For a long time, the computer experiment served as the best tool to investigate 
polymer systems. In this approach, the most widely used model was and is until now 
^ , the model of self-avoiding random walks (SAW) on a regular lattice. In spite of its 

simplicity, it perfectly captures the universal conformational properties of polymers in 
a solvent. Among the analytical descriptions of polymer systems the most successful 
are the field-theoretical approach [T], the direct polymer renormalization [3] and the 
traditional Flory theory [I]. The values of critical exponents and universal shape 
parameters, obtained within the frames of these theories, are in a good agreement 
with results of computer experiments. 
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Fig. 1. Schematic presentation of polymer chain in an environment with structural defects 
in form of parallel lines. 



The study of polymers in solvents in the presence of structural impurities plays 
a considerable role due to the importance to understand their behavior in colloidal 
solutions and near microporous membranes [7]. It also connects with the description 
of proteins inside the living cells which can be treated as disordered environments 
[5]. The behavior of polymers in disordered environments encountered controversies 
for a long time. During more than fifteen years there has been a wide discussion 
whether the presence of uncorrelated point-like defects leads to a change in universal 
properties |9ll()lllll2ll:^ll4linil(ill7ll8llf)l2()l2ll22) . This problem was first posed by 
Kremer showing that there exists a new value of the size exponent ([1]) only when 
the concentration of defects is above the percolation threshold. Later, this result was 
confirmed both analytically ^23^ and numerically |24I25I26I27I28I29I30] . 

The density fiuctuations of obstacles create large spatial inhomogeneities and pore 
spaces, which are often of fractal structure [3T]. A special case arises when fractal ex- 
tended defects of parallel orientation are present in a system, which may cause the 
anisotropy of environment (see fig. [1]). It is expected, that in this case there are 
two scaling exponents, governing the components of the size measure ([T]) in direc- 
tions parallel and perpendicular to the orientadion of the extended defects |32I33] . 
The influence of such disorder on the critical behavior of magnetic systems with Sd- 
dimensional defects of parallel orientation within frames of the spin m— vector model 
was analyzed in refs. [32i34i35j . It was shown that there are two characteristic corre- 
lation lengths in the system, one parallel and other perpendicular to the defects. As 
it is known [I] , the scaling properties of polymers in solution can be described within 
an m— vector model taking the de Gennes limit (polymer limit m 0).Unfortunatly, 
in the anisotropic case this approach does not give a satisfactory result. 

In the present paper, we aim to obtain a quantitative description of the scaling 
behavior of polymer macromolecules in environments with e^j-dimensional extended 
defects of parallel orientation (causing the anisotropy of environment), applying both 
numerical and analytical approaches. In the next Section, we study the special case 
of defects in form of randomly placed lines of parallel orientation (e^ = 1), applying 
computer simulations. In the Section 3, we develop an analytical description of the 
problem within the des Cloizeairx direct polymer renormalization scheme. We end 
giving conclusions and an outlook in Section 4. 
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2 Numerical analysis 

2.1 The Method 

We study the conformational properties of self-avoiding random walks on a regular 
lattice with extended impurities in the form of parallel lines (see fig. P), applying the 
pruned-enriched Rosenbluth method (PERM) [36 . The method combines the original 
Rosenbluth-Rosenbluth algorithm of growing chains [37^ and population control |38j . 
The n-th monomer is placed at a randomly chosen neighbor site of the previous 
{n — l)th monomer {n < N, where N is total length of the polymer). If this randomly 
chosen site is already visited by a chain trajectory or belongs to an impurity line, 
it is avoided without discarding the chain and the weight Wn given to each sample 
configuration at the n-th step is: 

n 

Wn^Y[mi, (2) 
1=1 

where mi is the number of free lattice sites to place the Ith monomer. 

The growth is stopped when the total length N of the chain is reached (or, at 
n < A^, if a the "dead end" without possibility to make the next step is reached), 
then the next chain is started to grow from the starting point that is chosen randomly 
every time. 

To derive appropriate properties of the chain we apply a two-step averaging: first 
average over configurations as usual and then averaging over different realizations of 
disorder. The configurational average for any quantity of interest then has the form: 

M M 

{{...)) = — J2w^i...), Z^^J^W^r', (3) 

^ k=l k=l 

where the summation is performed over the ensemble of all constructed A^-step SAWs 
(M ~ 10^ in our case). The disorder averaging is given in the form: 

ro-^E(-.-), (4) 

^k=l 

where p is number of replicas (we take p = 400 in our study) . 

The weight fluctuations of the growing chain are suppressed in PERM by prun- 
ing configurations with too small weights, and by enriching samples with copies of 
high-weight configurations. These copies are made while the chain is growing, and 
continue to grow independently of each other. Pruning and enrichment are performed 
by choosing thresholds W^f and , which are continuously updated as the simula- 
tion progresses. For updating the threshold values we apply similar rules as in |39l40j : 
= C(Z„/Zi)(c„/ci)^ and = 0.2W^, where c„ denotes the number of created 
chains having length n, and the parameter C controls the pruning-enrichment statis- 
tics; it is adjusted such that on average 10 chains of total length N are generated per 
each tour [ID] . 

2.2 Results and Discussions 

We perform simulations for chains length up to = 300 monomers to measure the 
end-to-end distances and up to TV = 100 monomers for the asphcricity value. 
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Fig. 2. End-to-end distance component perpendicular (left) and parallel (right) to the 
impurity lines as functions of chain length in double logarithmic scale at several values of 
defects concentrations. 



Let us introduce the position vectors r„ = {a;^, cc^, x^} of each n-th monomer 
(1 < n < N). In contrary to the the purely isotropic case ([1]), space anisotropy 
causes different scaling behavior for the end-to-end distance components in directions 
parallel (i?^ = {{x^ — xf)'^) and perpendicular = {{x]^ — xl^ + {xj^ ~ xfy) 

to the lines of defects (we assume, that the defects are extended in the x^ direction). 
Namely: 

"CRp - TV^.ii ^ J^^N^-^^ (5) 

with i^ii 7^ i^_L. 

Analyzing our results for the end-to-end distance components the parallel and 
perpendicular to the lines of defects (see Fig. 2), we can see that when concentration 
of impurity lines is small, there is a crossover between two types of behavior: one for 
short length ('^ 40 monomers) and the other for longer chains. We can interpret this by 
the fact that at small concentration of defects the short chain does not feel the presence 
of impurities (when averaged distances between impurity lines approximately equal 
the length of the chain). In such situations, to calculate the value of size exponents 
we have to take into account only the second part as it is closer to the asymptotical 
regime. 

As it was suggested, there are two different size exponents v\\ and iy±, which can 
be evaluated applying the least-square fitting of results presented in Fig. 2 to the form 



Concenfra(ion of impurities 



Fig. 3. Size exponents (01 at various concentration of defects. 
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Fig. 4. Aspheriity as function of the chain length at several different concentrations of 
defects (left). The change of aspheriity with increasing the defect concentration at fixed 
N = 100 (right). 



([S]). For example, in the case of concentration 20 % we have ly^ — 0.261 ± 0.005 and 
i^ll = 0.86 ± 0.07 which should be compared with corresponding pure value. Results 
for iy±, :/|| at different concentrations of impurity lines are given in Fig. 3. Whereas 
exponent gradually tends to the value of 1 as concentration of defects grows (the 
polymer chain extends in direction parallel to defects), the exponent vj^ is always 
smaller and gradually tends to zero. 

To analyze the shape of polymers we calculate the components of gyration tensor 
given by: 



N N 



Qij ~ y ^ y ^ (-^K ^m)(-^n -^m) (6) 



m— 1 n— 1 



where i,j — 1,2,3. The shape of a typical polymer chain conformation can be de- 
scribed by a quantity: 



Our results for {A3) are given in Figure 2] One can see, that at small concentrations 
of impurity lines the polymer chain does not change its shape in spite of the existence 
of two different characteristic lengths. Increasing the defect concentration to medium 
concentrations the shape become more prolate and tends to a rod-like state. 

As a result, polymers in an anisotropic environment with defects aligned in a given 
direction are elongated in this direction. In our case this is caused by the impurities 
that are extended along this direction and which repel the polymer. 



3 Analytical approach 

Let us start with a continuous model, where the polymer chain is presented as a path 
parametrized by r(s), with s — 0...S. An effective Hamiltonian of the system is 
given by: 

H=]^j^ (^) ' + " ds"5{r{s') ~ r(s")) ds + 



s 







V{v{s))ds. (8) 
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Note that r(s) is a d-dimentional vector with unit [length] while the parameters S, 
s have units [length]^. S is also called the Gaussian surface. Here, the first term 
describes the chain connectivity, the second term reflects the excluded volume effect 
with coupling constant u, and the last term arises due to the interaction between the 
monomers of the polymer chain and the structural defects in the environment given 
by potential V. 

We are interested in the case, when defects are correlated in Sd dimensions and 
randomly distributed in the remaining d—Sd- The correlation function of the densities 
of defects is thus assumed to be given by |33l32j : 



V(r(s))F(r(s')) = vS'^-'''{r{s) - r{s')). (9) 

Averaging the partition sum of the model over different realizations of disorder 
with ([9]) results in: 

Zdis = y X>{r} exp -ly^j^^^^ ds~uj\s' ds" 8{r{s') - v[s")) ds- 
[\s' r ds"S''-'^{rd^,,{s')-rd^,,{s"))ds , (10) 

Jo JQ 

Here, rd-t^{s) denotes the component of position vector r(s) perpendicular to the 
orientation of the extended defects and the functional integration J 'D{r} is to be 
performed over different chain configurations. Note that the coupling u must be pos- 
itive, which corresponds to an effective mutual repulsion of the monomers due to the 
excluded volume effect, whereas the coupling v should have negative sign, weakening 
the covolume effect. 

To evaluate quantitative estimates for the conformational characteristics of such 
a model, we apply the direct polymer renormalization scheme [3]. All properties of 
interest can be found in the form of a perturbation theory series in an coupling 
constants. In particular, for the components of the end-to-end distances ([5]) we found: 

, ^S{d~ed)\l + ^ + -, I, (11) 

d-.a y d)y (2-f)(3-f) (2-i^^)(3-^^)y ' ^ ' 

here Zu = u{2'k)-'^/'^S'^-'^/^ and z„ = ^;(27r)-(''-'^<*)/25.2-(d-e^)/2 ^^^.g dimentionless 
couplinds. Note, that £d indicates the dimentionality of the subspace parallel to the 
defects and d — Sd the dimentionality of the corresponding perpendicular subspace. 

In the asymptotic limit S oo the physical quantities, presented in the form 
of series expansions in the coupling constants, are, however, divergent. To obtain 
asymptotical values of the corresponding physical parameters a renormalization of 
the coupling constants needs to be performed '31 . The scaling exponents attain finite 
values when evaluated at the stable fixed point of the renormalization group trans- 
formation. The fixed points are defined as the common zeros of the fiow equations, 
which in our problem read: 

fiy ^ 5zy - %zl - 4:ZuZv, (13) 
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here e ~ A — d, 5 = e + Ed- We find four fixed points: 

< = 0, z:=0, (14) 

z:-e/8, < = 0, (15) 

z: = 0, z:=S/8, (16) 

z: = e/2-3<5/4, < - <5/2 - 6/4. (17) 



The first of them is the well known Gaussian fixed point, corresponding to an ideal- 
ized polymer chain (random walk) without excluded volume interactions. The second 
describes another well known problem namely the polymer in a good solvent (disorder 
is absent). The third point describes the Gaussian polymer in an anisotropic environ- 
ment and the last (most interesting one) describes SAW in anisotropic environment. 
Note, that at = we restore the situation with uncorrelated point-like defects, 
studied previously j23l41ll0ll7) . To correspond to a physical critical point of the sys- 
tem, the given fixed point should be stable and physically accessible. Unfortunately, 
in our problem the two last points, which are of main interest, are non-physical (e.g. 
attain appropriate values > and Zy < only in the unphysical region Sd < 0) 
and thus, cannot provide estimates of scaling exponents. The same problem appears 
for the coupling flow functions of the m— vector model with extended defects j35j . 
when one tries to take the m — > limit. 

A similar problem of the absence of stable and physically accessible fixed points 
exists also in the case of uncorrelated point-like impurities |10ll7j . but it was solved 
by adsorbing the interaction with disorder into the excluded volume interaction due 
to special symmetry [23]. However it does not work in our case of extended defects. 
In spite of the problems evaluating numerical estimates of the critical exponents, our 
analytical work confirms the existence of two characteristic lengths for polymer in 
anisotropic environments. 



4 Conclusions 

We analyzed the conformational properties of fiexible polymer chains in disordered 
anisotropic environments. The anisotropy is introduced in the form of e^-dimensional 
extended impurities of parallel orientation. Both numerical and analytical approaches 
are developed. 

In numerical part of our work, we studied the scaling properties of a SAW model 
on a lattice with defects in a form of lines {sd = 1), extended along a fixed coordinate 
direction. We conclude that there exist two critical exponents that describe the size 
measure of the polymer chain in directions parallel and perpendicular to the defects. 
The parallel exponent is always larger than the corresponding value for polymers 
in a pure solvent and gradually approaches the limit of 1 when increasing the de- 
fect concentration. The perpendicular exponent i'± is always smaller and reaches the 
limit of 0. We also found, that increasing the concentration of defects, the shape of 
the polymers becomes more anisotropic, elongated in the direction parallel to the 
extended defects and gradually reaches the limit of a rod-like shape. 

Our analytical study is developed on the basis of a continuous chain model with 
applying the direct polymer renormalization scheme. However, we encounter some 
controversies when analyzing stability and physical accessibility of fixed points, cor- 
responding to a critical point of the system. Clarification of this problem will be the 
subject of our further work. 
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